wo 2004/029568 PCT/AU2003/001280 

-1- 

TITLE 

Distributed temperature sensor using Optical Reflectometry 
INTRODUCTION 

5 

This invention relates to improvements in optical time domain 
reflectometry and more particularly relates to a distributed temperature sensor 
and its method of operation. 

10 BACKGROUND OF THF. TTWir.ivnrTOiv 

A distributed temperature sensor (DTS) measures temperatures 
along an optical fibre that is located in thermal contact with an object to be 
measured. The principle of operation of a DTS is analogous to that of RADAR and 

1 5 SONAR. In RADAR, the total time that an electromagnetic probing pulse takes to 
travel from the source to a distant object and that its reflection takes to return to 
the origin, together with the known speed of the electromagnetic wave, allows the 
location of a distant object to be deduced. In SONAR an acoustic probmg signal is 
employed. In DTS systems a very short pulse of laser light Oess than 100 ns) is 

20 used as the probe. After leaving the laser the Ught pulse travels through an optical 
fibre to an optical coupler and into the sensing fibre. As the light pulse travels 
along the sensing fibre, the pulse intensity is attenuated by scattering and 
absorption in the fibre material. 

25 Modem high-purity optical fibres have a wmdow of low absorption 

for wavelengths from approximately 0.8 to 1.7 ^m. Within this window, the optical 
losses are dominated by Rayleigh scattering, which contributes the following term 
to the total attenuation: 

30 where Xis wavelength, n is the refractive index, is the isothermal compressibility 
9tTfikis Boltzmann's constant and 3>is a constant related to the glass anneal 
temperature. The attenuation coefficient is usually expressed in units of dB/km. 
The probe pulse is also attenuated by mteractions with the optical medium that 
change the wavelength of the light. BrlUouin and Raman scattering are the most 

35 important of these relatively weak scattering mechanisms. These interactions 
cause light to be scattered back towards tihe proximal end of the fibre as light of 
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different wavelengths. The backscattered power received for Raman scattering at 
a particular wavelength K > X (known as the Stokes component) is given by 

while for \a < X(the anti-Stokes component) the backscattered power is given by 

^ ^""^ Xl[oxp{hv/kT)-l]' 

where h is Planck's constant, v is the frequency shift of the scattered light and Tis 
the temperature. Thus, while most light energy is transmitted in the forward 
direction along the fibre, a smaU fraction of it is scattered backwards, where it 
may be detected and analysed. Of the backscattered Ught, the Rayleigh scattering 
1 0 occurs at the probe wavelength and is relatively insensitive to temperature changes 
in the fibre, while the Raman scattering is shifted in wavelength from the probe 
light and has an explicit temperature dependence. 



In general, the accurate derivation of quantities based on the 
1 5 measurement of light intensity is most conveniently made through the 

measurement of light intensity ratios, because the absolute intensity (or power) is 
difficult to measure accurately. The intensity of an optical signal can be influenced 
by a large number of variables in addition to the quantity of interest For example, 
the power output of the source, the efficiency of the detector and the efficiency of 
20 the optical elements could all be affected by changes hi ambient temperature or 
humidity. Some components may be subject to ageing effects. These effects can be 
difficult to predict or identify and are therefore difficult to model reliably. 
Alternatively, the effects can be reduced through the provision of a stable thermal 
environment and suitable calibration means. In the art, it has been argued that 
25 various combinations of these approaches (ratiometric and calibration) can 

provide practical and efficient solutions to obtaining accurate measurements of 
temperature distribution, given all of the various sources of uncertainty that 
apply. 

30 The method of detection and analysis varies between different DTS 

embodiments based on glass optical fibres. In the earliest embodiments, a 
diffraction grating was used to filter out a band of backscattered wavelengths close 
to the laser wavelength (mainly the Rayleigh scattering). The Stokes and anti- 
Stokes Raman wavelengths were allowed to pass to separate detectors and the 
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intensity ratio of these components was used to derive the temperature as a 
function of range in the fibre (see GB 2,140^54A). 

An improved method was subsequently devised, whereby the 
5 Rayleigh scattering and anti-Stokes Raman scattering are selected for 

measurement by separate detectors (see GB 2,183,821A). These intensities are 
compared in a ratio device to give an indication of the temperatures in the fibre. It 
is claimed that this arrangement permits a much faster response than the prior art 
method, as the Rayleigh scattered light is much more intense than the Stokes 
1 0 Raman and can be sensed using relatively simple and inexpensive equipment 

In a further development, a method was devised whereby a single 
spectral band of the backscattered radiation (usuaUy a region of the broad anti- 
Stokes spectrum) is selected for analysis (see US 4,823,166). The method uses a 
calibration function to deduce the temperature distribution from the measured 
backscatter power. The data conversion may be carried out either using a 
tabulated variation of the backscatter factor with temperature, or via a theoretical 
model that relates absolute temperature to intensity as a function of \„. 



15 



20 



It is claimed that the embodiments described in US 4,823,166 remove 
the need for corrections to be made for the diflference in fibre attenuation between 
the Stokes and anti-Stokes wavelengths. It is also claimed that the system offers 
enhanced sensitivity to temperature changes, reduced sensitivity to drifts in the 
source wavelength and a simplified optical arrangement. Short-term changes in 
25 the energy and wavelength of the source can be detected and corrected by 
monitoring a short reference section of the fibre that is held at a constant 
temperature in a temperature-controlled chamber. However, since this approach 
relies on the accurate measurement of Intensity in a single spectral band. Its 
effectiveness is critically dependent on the elimination of variations m bJckscatter 
factor tiiat arise from non-temperature (NT) factors. In particular, axial variations 
in the fibre loss are of particular concern in the current context. These are 
typicaUy associated witii built-in or acquired defects in tiie fibre fliat cause 
temperature-independent variations fai flie scattering coefficient 



30 



35 



Specification US 4,823,166 suggests three ways in which the effects of 
variations in the fibre loss may be eliminated from the measured temperature 
distribution. The first metiiod involves performing tiie measurement from each 
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end of the optical fibre. The effects of any propagation losses are eliminated by 
calculating the geometric mean of the backscatter signals measured from both 
ends of the fibre and returning from a particular location. Unfortunately this 
approach adds to the instrumental complexity and is less convenient to deploy 
5 than a single-ended measurement arrangement 

A second approach involves calibrating the entire fibre before 
installation with a known temperature distribution. The sensor then measures 
departures of the backscatter intensity from those determined at the time of 
10 calibration and hiterprets them in terms of a temperature variation. However, this 
approach restricts the system to use with fibres for which a calibration has been 
performed and requires recalibration if the fibre properties change. The third 
approach makes provision for the removal of the filter to faciUtate measurement of 
the total backscatter signal in the reference section, or over the entire fibre length, 
15 so that a normalisation can be performed. The total backscatter signal is 

dominated by the Rayleigh scattering, which is relatively temperature insensitive, 
but sensitively reflects the fibre loss characteristics. However, the need to remove 
the filter to perform the normalisation procedure adds to the complexity of the 
optical system and remains a drawback for practical operation. 



20 



It is these issues that have brought about the present invention. 
SUMMARY OF Tmr. INVENTION 



I" accordance with one aspect of the present invention there is 
provided a temperature sensing method In which pulses of optical radiation are 
launched by a laser diode into an optical fibre and optical radiation backscattered 
from the fibre is detected, the method comprising passing the backscattered 
radiation through a shigle optical filter whereby a first signal is recorded at the 
30 anti-Stokes Raman wavelength from a signal launched by tiie laser diode In a laser 
mode and a second signal is recorded at the Rayleigh wavelengtii from a signal 
launched by the laser diode In a Hght emitting diode mode, and a comparison is 
made of the two signals to provide an indication of temperature. 



35 



Preferably, the comparison produces the quotient of the anti-Stokes 
Raman wavelength divided by the Rayleigh wavelength. 
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Preferably, a photomultiplier tube is used to count photons through 
the single optical fUter, the photon count providing an indication of temperature 
and the arrival time of the signals allowing identification of the position of that 
temperature along the length of the optical fibre cable. 

5 

In a preferred embodiment the laser diode, the single optical filter 
and the photomultiplier tube are coupled to the optical fibre via connectors and at 
least one temperature sensor is positioned downstream of the connectors to correct 
error signals emanating from flie connectors. 



10 



The method further comprises stabilising and controlling the 
temperature of the laser diode and photomultiplier tube. 



Preferabfy, a short reference section of flie optical fibre is held at 
1 5 constant temperature in a temperature-controlled chamber. 

In a preferred embodiment, two distinct controUed temperature 
regions are provided and the temperatures are monitored through the use of 
sensors in each region, and the method monitors the efifect of the temperature- 
20 controlled region on the accumulated photon count 



Preferably, the pulses of optical radiation are transmitted at a 
wavelength between the ultra-violet and infrared spectrum and particularly at 
wavelength in tiie range 775-800nm. 

25 

In a preferred embodiment, the laser diode operates at a 
transmission power of less than IW. 



In accordance with a fiirther aspect of the present invention there is 
provided a distributed temperature sensor comprising a laser diode adapted to 
launch pulses of optical radiation into an optical fibre arranged to be located in 
thermal contact with an object, and a single optical filter to detect optical radiation 
backscattered from the fibre, the distributed temperature sensor being adapted to 
operate in accordance with the method described above. 



Preferably, the single optical filter is coupled to a photomultiplier 

tube. 
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BRIEF DESCRIPTION OF THE DRAWTNrtS 



An embodiment of the present invention will now be described by 
5 way of example only with reference to the accompanying drawings m which: 
Figure 1 is a schematic illustration of the basic components of a 
distributed temperature sensor. 

Figure 2 is a graph of light spectrum against wavelength illustrating 
the major spectral features and location of an optical filter, 
^® Figure 3 is a graph iUnstrating an accumulated photon count, 

Figure 4 is a schematic illustration of a distributed temperature 
sensor and temperature calibration means. 

Figure 5 is an illustration of an accumulated photon count usmg the 
calibration means of Figure 4, 

^ ^ Figure 6 is a circuit dii^ram relating to the LD power supply. 

Figure 7 is a circuit diagram relating to the LD driver electronics, 

and 

Figure 8 is an illustration of the operational amplifier (Op-Amp) 
output of the LD driver circuit with laser diode voltage and current waveform 
20 signals of approximately 10 ns width. 

In the distributed temperature sensor (DTS) 10 shown in Figure 1 a 
very short pulse of laser light (less than 100 ns) constitutes a probe. After leaving 
the laser 11 the probe light pulse travels through an optical fibre 12 to an optical 

25 coupler 13 and into the sensing fibre 14. As the probe pulse travels along the 
senshig fibre, the glass atoms absorb or scatter a smaU fraction of the incident 
light and the pulse intensify is attenuated. The fraction of the hicident light that is 
scattered back along the fibre wiU exhibit a spectrum such as that shown in Figure 
2. In the DTS»s normal operating mode the bandpass filter 15 aUows only the 

30 shorter wavelength anti-Stokes light to be detected. The bandpass filter transmits 
a selected wavelength range, as shown in Figure 2, to the photomultiplier light 
detector 16. 



The DTS operates in a photon counting mode. The intensity of the 
light returning to the detector is sufficiently low as to comprise mostly of separate 
photons arriving at distmguishable times. The detection system represents each 
photon as an electrical signal that is digitaUy recorded as an event labeUed with the 
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10 



35 



time of arrival of the photon, measured from when the laser pulse was fired. Light 
to be recorded is selected with a single filter set at a central wavelength of 750-755 
nm (for a probe wavelength of 780-785 nm) and a full width half maximum pass 
band of 10-20 nm, a transmission of greater than 75% in the passband and a 
rejection of better than lo". The digital nature of photon countmg has the 
advantage of reducing the sensitivity of the detection system to environmental 
variations. On the other hand, photon counting allows sensitive detection of 
extremely low fight levels, which aUows the use of a relatively low power laser 



source. 



The total time taken for fight to travel from the source to the point of 
scatter and for the scattered fight to travel back to the detector are recorded at 
discrete intervals of 2 ns (bins). It is the intensity of the returned fight together 
with its time of arrival that aUows the system to deduce the temperature profile 

15 along the sensing fibre. The process is repeated many thousands of times per 
second and the results are accumulated over a predefined period (lypicafiy 
minutes). Figure 3 shows a typical accumulated photon count result The 
accumulated photon count as a function of position along the fibre is then 
downloaded to a central processing unit where the temperature distribution is 

20 calculated. The DTS requires high-speed electronics for the LD driver and photon 
detection circuits but uses the essential advantages of the photon counting method 
to simplify aU other aspects of DTS operation. 

Although the use of photon counting and a relatively low power laser 
25 are inherenfly favourable for reducing tiie sensitivity of the system to 

environmental and qflier variables, it remains necessary to stabfiise and control 
flie temperature of the laser diode and photomultipfier tube. This serves to reduce 
the sensitivity of the system to drifts in tiie source wavelengtii with temperature 
and reduces the background noise level, respectively. Further aUowance for short- 
30 term changes in tiie energy and wavelengtii of the source can be made by 
monitoring a short reference section of the fibre fliat is held at a constant 
temperature in a temperature-controlled chamber. One means of arrangtag such a 
reference section is filustrated in Figure 4, where temperature sensors Trl and Tr2 
monitor two distinct temperature regions. The effect of this temperatiire 
caUbration region on tiie accumulated photon count is shown in Figure 5. 

In summary the DTS monitors the backscattered signal using: 
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• the photon counting method to record light intensify, and 

• the arrival time to identify position. 

5 The DTS distributed sensor 10 is characterized by: 

1) A laser diode light source 11. 

2) A photomultiplier light detector 1 6. 

3) A 50/50 optical coupler 13. 
10 4) An optical filter 15. 

5) An optical fibre 14 in which temperature is measured. 

6) Temperature calibration system. 

7) Analogue control circuitry. 

8) High-speed laser driving circuitry. 

15 9) High-speed photon detection circuitry. 

10) High-speed photon counting and accumulation circuitry. 

11) A CPU for control, data collection and processing. 

The Ught probe 11 is generated by a Sanyo DL7140 785nm, SOmW 
20 continuous wave output, single-mode, AlGalnP, index guided structure, solid-state 
laser diode. The electronic drive circuitry is shown in Figure 6. The drive circuitry 
aUows for two LD operating modes: LASER and LEDing mode. A complete 
description of these modes is given below. 

^5 Figure 2 identifies three spectral components: the shorter 

wavelengtii anti-Stokes Raman (AS), Rayleigh (R) and the longer wavelength 
Stokes Raman (S). All these spectral components contain to varying degrees both 
temperature (T) and non-temperature (NT) information. The DTS challenge is to 
recover only tiie temperature information. Conventional distributed temperature 

30 tiieory provides that a temperature determination that is essentiaUy independent 
of the KT effects can be made by taking the ratio AS/S or alternatively the ratio 
AS/R. The ratio of the amounts of light backscattered into the two Raman bands is 
given by: 

_ Anti - Stokes Intensity ( ft^^ 

Stokes Intensity ~ /Ll,^^\~'kT) 

35 and is tiius a function of the temperature of the fibre at the scattering site. 
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In such cases measurements of light at two separate wavelengths are 
required. In some embodiments of this kind of technology aU three spectral 
components are detected and recorded. In these cases there is an increase in 
5 system complexity and more components are required, leading to increased cost. 
In addition, complicated corrections have to be made for the fact that the 
absorption loss in the fibre differs slightly at the two wavelengths and is also 
temperature dependent 

The DTS system described hereunder aUows certain non- 
temperature dependent effects on the optical fibre to be measured in order to 
remove these effects from the final result. The NT effects are measured by means 
of the Rayleigh signal, which is relatively temperature insensitive. This procedure 
exploits the characteristics of a laser diode whereby at low current the device 
behaves as a light emitting diode whilst at higher currents it has the properties of a 
laser. 



15 



The electroluminescence of the semiconductor laser and the light- 
emitting diode LED result from a current flow through a p-n junction to which a 
20 voltage has been applied. Recombination of the carriers injected across the 

junction results in the emission of Ught A measure of the quality of this process is 
the quantum efBciency, defined as the ratio of the number of emitted photons to 
the number of electrons crossing the p-n junction. The spontaneous emission of 
LEDs is characterised by low quantum efBciency and a relatively broad spectrum. 



25 



In the semiconductor laser both the electrons and the radiative 
emissions are confined by a cavity with partially reflective surfaces to produce 
stimulated emission. Stimulated emission in such an arrangement results when a 
photon with an energy slightly greater than the energy gap can interact with an 
30 electron in the conduction band and cause the electron to recombine with a hole in 
the valence band. This recombination process results in the emission of a photon 
identical to the photon that caused the recombination process, and the number of 
photons is mcreased. At low current, the laser light intensity is smaU and results 
from spontaneous emission, as in a LED. However, when the stimulated emission 
exceeds the internal losses, the laser threshold is reached and the light output rises 
rapidly with the current. Above the threshold current, most of the current flowing 
into the p-n junction results in laser emission, the quantum efficiency is much 



35 
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higher than for an LED and the emitted light is almost, but not quite 
monochromatic This dual nature of semiconductor lasers is utilised in the present 
embodiment. 



5 The principle of operation is: 



10 



15 



20 



25 4. 



a) Make a measurement with LD in LASER mode yielding a record of 
temperature and non-temperature data. 

b) Make a measurement with LD in LED mode. In this case light 
passing through the filter is predominantly due to the Rayleigh effect 
and is not temperature sensitive. 

c) MathematicaUy deduce a record from both signals that only reflects 
the temperature dependence. 

This functionality is achieved through a novel circuit, consisting of 
regulated power supply and laser driver elements, shown in Figures 6 and 7 
respectively. 

The features of the circuit are: 

!• Operation in LASER mode. 

2. Operation in LED mode. 

^- The probe pulse can be widened to increase the power so that better 

temperature resolution can be obtained. 
An ability to control the current so as to produce a high-speed 
probing pulse for better range resolution. 



Power Supply: 

^® When the equipment is first turned on there is no power to the LD 

drive circuit A 5V supply (labelled «+5V slow" in Figure 6) is appUed to the LD 
drive circuit a few seconds after other system circuitry is stable. The voltages Vh 
and -V are variable and are used to control aspects of the laser operation as 
described below. 

35 



Quiescent circuit condition: 
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In the quiescent state, a negative bias voltage (-V) is applied to the 
LD to provide a level of current just below the value required for light emission. 
Its function is to reduce the amount of current required to create light emissions 
and thereby reduce the magnitude of the voltage swing and slew rate required of 
5 the operational amplifier (Op-Amp). 



The LD is biased, in its normally-ofiT state, by being connected 
between the Op-Amp output (normally OV) and -V. The bias voltage -V is chosen 
to be about -1.6 V to maintain the LD just below its turn-on point It is critical to 
1 0 keep the current very low in this state as modem heterostmcture type LDs may 
act as reasonably efBdent LEDs even at sub-milUamp levels. Further bias 
adjustment is provided by a resistor from +5V to the hiverting input, particularly 
relevant for the period between steps 3 and 4 below. 



15 Dynamic behaviour: 



1 Dynamic events commence when a nominally 50 ns wide TTL positive- 
gomg trigger pulse is applied to the circuit. The shape of this pulse Is not 
particularly critical and it does not define the circuit's output pulse width. 
20 This trigger takes two paths. 

a. The trigger signal is reshaped and delayed. The input trigger pulse is 
reshaped into pulses with faster rise times by the 10 Schmitt triggers 
connected in series and internal to the 50A-10250 TTL digital delay 
line integrated circuit package. Reshaped versions of the input 
trigger signal appear at successive outputs of this package, each of 
which is delayed by 2.5 ns from the previous one. 

b. To improve the LD»s tum-on-time a small bias current is applied to 
the diode prior to the application of the trigger signal from the delay 
line. This bias current is due to a small fraction of the input trigger 

30 P«lse that is appUed to the non-inverting input of the operational 

amplifier (HFA1130). 

2. In the circuit shown in Figure 7, one of the reshaped signals delayed by 2.5 
ns (Tl) is connected to the non-Inverting input of an HFA1130 Op-Amp. At 
35 this point the inverting input source is still zero, and the amplifier output is 

large and positive. Current limit into the LD is provided by a combination 
of the Op-Amp Vh output limit and a 15 O series resistor. 
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It should be noted that the HF1130 Op-Amp has two critical differences 
from other so-caUed "fast" Op-Amps. Firstly the HF1130 contains high and 
low output voltage limit circuitry (Vh and VL inputs), which prevent the 
5 internal circuitry from saturating. This allows very large gains to be used, 

which effectively improve the rise time of the output without the normaUy 
associated dampened response caused by saturation in other Op-Amps. 
Secondly, the HFA1130 output drive contains special high speed/current 
components with more in common to a Norton Amplifier than an Op-Amp. 
10 The output rise time is 6000 VZ/ts which is much faster than typical driver 

circuits with speeds of only 100-500 Wlpis for low loads and small signals. 
All Op-Amp inputs comply with manufacturers recommendations for input 
impedance and circuit layout as far as practicable. 

15 3. The second delayed trigger pulse from the delay line is connected to the 
inverting input of the Op-Amp. Any one of the delayed outputs (for 
example T7) can be used. The final output pulse width is determined by the 
time difference (for example T1-T7). This signal causes an inverted and 
greatly amplified output which reverse biases the LD and completety 

20 switches it off much faster then zero biasing would. The Op-Amp inputs 

have a common mode voltage of approx 2.5 V at this stage. 

4. When the first delay line output returns to zero volts, but before the second 
output does, the diode (U16) between the delay line outputs conducts and 
25 briefly shorts them. This begins to switch off the LD reverse bias and 

removes an output spike that would otherwise be present. A smaU signal 
silicon diode, such as the common IN914B, with high speed under these 
conditions is most appropriate. 

30 5. Both delay line outputs have returned to zero volts and LD is returned to 
non-conducting state. 

Since the voltage applied to the LD resistor circuit does not directly 
indicate the magnitude of the laser current, the LD current pulse is deduced from 
35 the potential across the resistor. This method better predicts the optical power 
output, the turn-on rise time, and makes visible the stored charge flowing in 
reverse direction from the LD during turn-off. Typical behaviour is illustrated in 
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Figure 8. 

Operation in LED mode: 

5 In order to switch the output of the LD to the broad LED spectrum, 

the power supply circuitry shown in Figure 6 includes a software-controUed switch 
that connects Vh to a variable resistance VR3. This resistance acts as a voltage 
divider and is used to reduce Vh to a level of about 250 mV, thereby reducing the 
current into the LD. The LD is now unable to reach the laser threshold and the 
10 light emission Is predominantly spontaneous, as in an LED. 

In all other respects, the dynamic behaviour of the LD in LED mode 
is the same as that described above. The broad LED spectrum now includes a 
component that coincides with the filter window, aUowing the Rayleigh scattering 
at the filter wavelength to be detected directly, using the same optical arrangement 
as that used to measure the Raman scattered light The reduced intensity of light 
emitted at the filter wavelength by the LD in LED mode is compensated to some 
extent by the higher efficiency of the Rayleigh scattering, relative to the Raman 
scattering. 

20 

In this way, the temperature-dependent Raman scattering and 
temperature-independent Rayleigh scattering can be obtained using a common 
optical arrangement. AnomaUes in the fibre transmission (associated for example 
with kinks or strain In the fibre) can therefore be identified qualitatively in the 
Rayleigh scattering. These anomaUes may then be discounted in the temperature 
dependent Raman measurement Alternatively, the Rayleigh scattering 
hiformation can be used as the basis for a mathematical correction of the 
temperature distribution obtained from the Raman measurement The correction 
must take into account the difference in attenuation for probe light at the LD 
30 wavelengfli and at tiie filter wavelength. This difference in attenuation is readily 
obtained by comparison of Rayleigh and Raman measurements taken with the 
fibre subjected to a known temperature distribution. 



15 
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35 



It wiU be appreciated that a number of alternative means for 
generating LED output from a LD wiU be apparent to those skilled in the art 



